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ABSTRACT 



Context. Studying the density profiles of galaxy groups offers an important insight into the formation and evolution of the structures 
in the universe, since galaxy groups bridge the gap between single galaxies and massive clusters. 

Aims. We aim to probe the mass of SL2S J02140-0535, a galaxy group at z = 0.44 from the Strong Lensing Legacy Survey (SL2S), 
which has uncovered this new population of group-scale strong lenses. 

Methods. We combine strong lensing modeling and dynamical constraints. The strong lensing analysis is based on multi-band 
HST/ACS observations exhibiting strong lensing features that we have followed-up spectroscopically with VLT/FORS2. To con- 
strain the scale radius of an NFW mass profile that cannot be constrained by strong lensing, we propose a new method by taking 
advantage of the large-scale dynamical information provided by VLT/FORS2 and KECK/LRIS spectroscopy of group members. 
Results. In constrast to other authors, we show that the observed lensing features in SL2S J02140-0535 belong to different background 
sources: one atz= 1.7 ± 0.1 (photometric redshift) produces three images, while the other at z = 1.023 ± 0.001 (spectroscopic redshift) 
has only a single image. Our unimodal NFW mass model reproduces these images very well. It is characterized by a concentration 
parameter c 2 oo = 6.0 ± 0.6, which is slightly greater than the value expected from ACDM simulations for a mass of M 2 oo ~ 1 x 
10 14 Mq. The spectroscopic analysis of group members also reveals a unimodal structure that exhibits no evidence of merging. The 
position angle of the halo is 8 = 111.6 ± 0.2, which agrees with the direction defined by the luminosity contours. We compare our 
dynamic mass estimate with an independent weak-lensing based mass estimate finding that both are consistent. 
Conclusions. Our combined lensing and dynamical analysis of SL2S 102140-0535 demonstrates the importance of spectroscopic in- 
formation in reliably identifying the lensing features. Our findings argue that the system is a relaxed, massive galaxy group where 
mass is traced by light. This work shows a potentially useful method for constraining large-scale properties inaccessible to strong 
lensing, such as the scale radius of the NFW profile. 

Key words. Gravitational lensing: strong lensing - Galaxies: groups - Galaxies: groups: individual (SL2S J02140-0535) 



1. Introduction 

In the context of structure formation and evolution, galaxy 
groups play a key role, spanning the regime between individual 
galaxies and galaxy clusters. Therefore, detailed studies of this 
intermediate regime of the mass spectrum are relevant. Galaxy 
groups have been mostly studied using X- ray and optical trac- 
ers, a s well a s numerical simulations (e.g . iHelsd on & P onmanl 
200f j I20031 lOsmond & Ponmanl l2004t I Willis et all I2005 



Finoauenov et alJ | 2007t iRa smussen & Ponma 



2007 



Mamon 2007; Faltenbacher & Mathews 2007; Yang et al. 
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* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de l'Univers of the 
center National de la Recherche Scientifique (CNRS) of France, and the 
University of Hawaii. This work is based in part on data products pro- 
duced at TERAPIX and the Canadian Astronomy Data center as part 
of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative 
project of NRC and CNRS. Also based on HST data as well as Keck 
(LRIS) and VLT (FORS 2) data. 



2008; van den Bosch e t~all l2008t fG astaldel lo et all l2007b 
D'Onghia et alJ I2005I; ISommer-Larsenl l2006i and references 
therein). Only recently, thanks to the unsurpassed combined 
depth, area, and image quality of the Canada France Hawaii 
Telescope Legacy Survey (CFHTLSfl a new population of 
lenses has been uncovered, with Einstein radii between ~ 3" 
and ~ 8", generated by galaxy-group-scale dark-matter haloes 
(this subsample of the CFHTLS forms the Strong Lensing 
Legacy SurvejQ (SL2S). ICabanac et al.ll2007l) . This new pop- 
ulation effectively bridges the gap between single galaxies 
and massive clusters, opening a new window of exploration 
in the mass spectrum. The first representative sample of these 
lenses was presented by Limousin et al. (2009). Since then, 
a detailed investigati on of some gro ups has been carried out: 
SL2SJ02 176-05 13 (lTuetal.1 120091) and SL2S J08544-0121 
(iLimousin et alJ 12010b . In this paper, we present a detailed 
analysis of SL2S J02140-0535 at z = 0.44, combining the 
dynamics of the group members, and both the strong and weak 
gravitational lensing. 



1 http : //www . cfht . hawaii . edu/Science/CFHLS/ 

2 http://www-sl2s.iap.fr/ 
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Fig. 1. The group SL2SJ02 140-0535 at z spec = 0.44. Upper left: composite CFHTLS g, r, i color image (1 square arcminute = 
341 x 341 kpc 2 ). Upper right: composite HST/ACS F814, F606, F475 color image (22" x 22" = 125 x 125 kpc 2 ). Note how the 
space-based image significantly improves the identification of multiple images and heightens the resolution allowing us to detect the 
substructure. Lower: CFHTLS i band (7.2' x 7.2' = 2.46 x 2.46 Mpc 2 ) centered on the strong lensing system. We draw in green the 
luminosity contours (corrected for passive evolution) corresponding to 10 6 , 3 x 10 6 , and 10 7 L.Qkpc~ 2 , respectively. Red rectangles 
show the location of the slits used for spectroscopy for the 16 confirmed members of the group and the magenta circle denotes a 1 
Mpc distance from the center of the group. 
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Table 1. Photometric and spectroscopic data of the galaxies and the arcs. 



ID Name 
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u* - g' 


g> -r> 


r' — i' 


/' - z? 
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0o 
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G\ 


0.00 


0.00 


2.35 ±0.11 


1.61 ±0.02 


0.86 ± 0.01 


0.36 ±0.01 


0.158 


-41.0 


0.4446 ± 0.0002 


G2 


1.16 


3.09 


2.99 ± 0.23 


1.64 ±0.03 


0.87 ± 0.01 


0.37 ± 0.01 


0.424 


64.2 


0.4449 ± 0.0003 


G3 


1.37 


1.25 


1.81 ±0.15 


1.47 ±0.06 


0.88 ± 0.04 


0.33 ± 0.03 


0.177 


18.26 


— 


A 


— 


— 


0.73 ± 0.09 


0.30 ± 0.09 


0.47 ±0.10 


0.31 ±0.11 


— 


— 


[1.7 ±0.1] 


B 


— 


— 


0.60 ±0.11 


0.07 ±0.11 


0.39 ±0.14 


0.36 ±0.14 


— 


— 


[1.7 ±0.1] 


C 


— 


— 


0.18 ±0.09 


0.58 ± 0.09 


0.60 ± 0.09 


0.41 ± 0.09 


— 


— 


1.023 ±0.001 



Notes. Column (1) lists the identification for each object as in Fig.Q] Columns (2) and (3) list the relative positions measured with respect to the 
brightest group galaxy (BGG) with coordinates a = 02:14:08:097 and 6 = -05:35:33.77. Columns (4-7) list the CFHTLS colors. Columns (8) and 
(9) show the geometric parameters: ellipticity e and position angle 9 , which gives the orientation of the semi-major axis from the horizontal line 
in the image, measured counterclockwise. In the last column we list the redshift; values in square brackets are photometric redshifts. 



Some researchers have demonstrated the advantages of com- 
bining the different lensing regimes with the constraints that 
come from the stellar kinem atics of the centra l brigh test cluster 
galaxy (BCG). For example. [Sand et all d2002ll2004 combined 
strong gravitational lensing with the velocity dispersion profile 
of the BCG and found th at the model s favored a shallower slopes 
in the NFW profiles. Newm an et alj (|2009) following the above 
works and adding weak-lensing analysis, sampled the dark mat- 
ter profile over three decades in radius in the lensing cluster 
Abell611. Despite all these works, on group scales the efforts 
in that direction are just beginning. First, because strong lens- 
ing in groups is a relatively new field of research, and second, 
owing to observational limitations, i.e., there are fewer lensing 
arcs in groups and detecting weak lensing by groups of galax- 
ies is challenging. McKean et al. (2010) probed the inner part of 
a group-scale h alo using strong gravita tional lensing and stellar 
kinematics, and lThanjavur et all ([2010) showed that is possible 
to combine strong lensing and galaxy dynamics to character- 
ize the mass distribution and the mass-to-light ratios of galaxy 
groups. 

SL2S J02140-0535 was previously studied bv I Alardl d2009l) 
using a perturbative reconstruction of the gravitational lens. He 
found that the local shape of the potential and density of the 
lens, inferred from the perturbative solution, reveals the exis- 
tence of an independent dark component that does not follow 
light. He argued that this particular shape of the dark halo is due 
to the merging of cold dark matter halos. Given these findings, 
the group properties represent a challenge, and to describe its 
mass distribution and perform a detailed study of the object is 
of significant worth and deserves attention. Our first aim in this 
work is to build a strong lensing mass model for SL2S J02140- 
0535 using lensed features observed in HST/ACS data combined 
with ground-based spectroscopy. The second is to combine the 
strong lensing constraints with the dynamics of the groups mem- 
bers to study the mass profile out to the scale radius. Finally, we 
wish to compare the masses obtained with the three methods, 
namely, strong lensing, weak lensing, and group dynamics, to 
shed some light on SL2S J02140-0535. 

The paper is organized as follows. In Sect.|2] we describe 
the observational data images and spectroscopy and present the 
multiple image identification. In Sect.[3] we define how the dy- 
namical constraints (on galaxy and group scales) enter into our 
strong lensing model. We also depict the profiles for both, the 
group and the galaxy-scale mass components. In Sect.[4] we 



present the strong lensing as well as the weak lensing analy- 
sis of SL2SJ02 140-0535 and the main results of the work. In 
Sect.|5] we summarize and discuss our results. Finally in Sect.|6] 
we present the conclusions. All our results are scaled to a flat, 
ACDM cosmology with Om = 0.3, Qa = 0.7 and a Hubble con- 
stant Ho = 70 kms -1 Mpc~'. All images are aligned with WCS 
coordinates, i.e. north is up, east is left. Magnitudes are given in 
the AB system. 



2. Observations and multiple images identification 

SL2S J02 140-0535 has been imaged by both ground and space 
based telescopes. From the ground, this group was observed 
in five filters (u*,g',r',i',z') as part of the CFHTLS using 
the wide field imager MegaPrime, which covers ~ 1 square 
degree on the sky, and a pixel size of 0.186". From space, 
the lens was followed-up with the Hubble Space Telescope 
(HST). Observations were performed in snapshot mode (C 15, 
P.I. Kneib) in three bands with the ACS camera (F814, F606, 
and F475). Moreover, a spectroscopic follow-up of the arcs and 
the group members have been carried out from VLT and Keck. 



2.1. Imaging 

The strong lensing deflector in this group is populated by three 
galaxies (see Fig.[TJtop left). The lensed images consist of three 
arcs surrounding the deflectors: one arc (labeled A) situated 
north of the deflector composed by two merging images, a sec- 
ond arc in the east (labeled B), and a third one (labeled C) in 
the south. In this ground-based image, we see a possible system 
consisting of ABC images that seem to have compatible colors. 
In Fig.[TJ(top right), we show a color composition from the HST- 
ACS camera for the central part of the group. Since in the HST 
image, arc C exhibits different color from both A and B, is it pos- 
sible that this arc is not associated with the proposed multiple- 
image system? This scenario is supported by the substructure 
in arcs A and B th at are not pr esent in arc C (see top right of 
Fig.Q]). Given that Alardl (I2009I) published an analysis based on 
this southern image being part of the multiply imaged system, 
special care will be used to properly assess the multiple image 
identification, which is crucial for the analysis. 

The photometry for the central galaxies labeled Gl, G2, and 
G3 (see top right of Fig. Q}, and for the arcs was performed in all 
CFHT bands with the IRAF package apphot. For the galaxies, 
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the magnitudes were derived by measuring fluxes inside a fixed 
circular aperture of 12 pixels (~ 2.22") for galaxies Gl and G2, 
and 6 pixels (~ 1.12 ") in galaxy G3. For the arcs, we employed 
polygonal apertures to obtain more accurate measurements. The 
vertices of the polygons for each arc were determined using the 
IRAF task polymark, and the magnitudes inside these apertures 
were calculated using the IRAF task polyphot. 

The general properties of the central galaxies as well as the 
arcs are presented in Table Q]. Column (1) lists the identification 
for each object. Columns (2) and (3) list the relative positions 
measured with respect to the brightest group galaxy (BGG) with 
coordinates a = 02:14:08:097 and 6 = -05:35:33.77. Columns 
(4-7) list the colors u* - g', g' — r', r' —i' and i' -z' ■ Columns (8) 
and (9) show the g eometric parameters derived from SExtractor 
(iBertin & Arno uts 1996), i.e., ellipticity e and position angle Go, 
which imply the orientation of the semi-major axis from the hor- 
izontal line in the image, measured counterclockwise. Finally, 
column (10) lists the redshifts (see the next subsections). We 
note in Table [TJthat the colors u* - g' and g' — r' for arc C dif- 
fer significantly from the colors in A and B, which have almost 
similar values (within the errors). 



2.2. Spectroscopy 

Two spectroscopic runs have bee n performed in t he field of 
SL2SJ02140-0535. Using LRIS (lOke et all [1991 on Keck 
(dichroic 680 and 560 nm, exposure 600s and 1200s with a slit 
of 1.5"), we measured the redshifts of one of the brightest galax- 
ies populating the strong lensing deflector (G2) and we found z 
= 0.4449 + 0.0003. 

We also used FORS 2 on VLT (P.I. V. Motta) with a medium 
resolution grism (GRIS 600RI) to target both group members 
and the strongly lensed features. Two masks were used and two 
exposures of 1400s were made with each mask. We selected the 
targets (other than strongly lensed features) as a two step pro- 
cess. First we choose galaxies whose g - i color is within 0. 15 of 
the g - i color of the t wo brightest g alaxies populating the strong 
lensing deflector (see. Limousin et al. 2009). Since the field was 
too crowded in the center but not fully sampled at larger dis- 
tances, we also selected all bright galaxies in the outer field re- 
gardless of their color. Fig.Q~]shows the slit configuration for all 
confirmed group members. Full details of the spectroscopy of 
the group members will be presented in a forthcoming publica- 
tion (Munoz et al., in prep.). Here we concentrate only on the 
spectroscopy of the lensed features. 

No spectroscopic features were found in arcs A and B. On 
the other hand, the spectra of feature C reveals a strong emission 
line at 7538.4 A. As we can appreciate in Fig. [2] the emission line 
has a spatial extension that covers the whole arc (although the 
slits do not completely cover arcs A and B, they cover the bright- 
est parts, and in arc B only a small fraction is outside). From this 
analysis, two possible scenarios emerge: system AB and arc C 
come from two different sources, or belong to the same source 
but the fluxes in system AB are not strong enough to detect the 
line. From the magnitudes in u* band, it is possible to measure 
the flux ratios for the arcs, fs/fc ~ 0.5 and /a//c ~ 2.5 (we used 
half of the flux in arc A to make the calculations). Assuming that 
the magnification is not very different in the arcs and they come 
from one single source, if the emission line is detected in arc C 
then it should be detected in arc A. It should even be detected in 
arc B, since the flux is only half of the flux in C and the feature 
is strong enough (see top and bottom panel of Fig. |2). That this 
is not detected, excludes the second scenario and supports the 



7600 




Fig. 2. Top panel. Spectrum of arc C. It shows a strong feature 
probably associated with [OII]/13727. Middle panel. VLT pre- 
image of the central part of SL2S J02 140-0535 with the slits po- 
sitions. Bottom panel. The 2D spectra of arcs A, B, and C in the 
same wavelength region. Note the extended emission line in C 
spectrum. 



conclusion that system AB and arc C are produced by different 
sources. 

The line detected in C may correspond to either [OII]/13727 
at z ~ 1 or Lya at z ~ 5. However, the low redshift solution 
is much more likely, and as we discuss below, the photometric 
redshift estimate supports this conclusion. Thus, assuming emis- 
sion from [OII]/l3727 and applying a Gaussian fitting, we obtain 
Zspec - 1-023 ± 0.001. This line coming from the whole arc C 
is incons i stent w ith the model of SL2S J02140-0535 presented 
in lAlardl (J2009). He assumed that arcs A, B, and C belong to 
the same source and he argued that the differences in colors ob- 
served in the HST-ACS images were related to the crossing of 
the caustic line by part of the source. His model predicts that 
one half of image C comes from the same part of the source as 
images A and B, but the other half corresponds to an area of the 
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Fig. 3. Output photometric redshift PDF (see text) for arcs la- 
beled A, B, and C (top to botton). The dashed vertical line corre- 
sponds to the spectroscopic redshift inferred assuming the emis- 
sion line comes from [OII]/13727. 



source that has no counterpart in the large arcs (A and B). The 
lack of spectroscopic features in arcs A and B does not support 
this scenario. 

As a matter of fact, it is known that arcs produced by one 
single source could have different colors if the source is located 
over a caustic line. For example, in the case of a spiral galaxy for 
which the disk (but not the bulge) is located over a caustic line, 
part of the disk can be mapped onto an arc with colors that do 
not agree with that of the counter-image. In some extreme cases, 
a tangential arc can indeed show no evidence of having been 
formed by two mirror images (e.g. Smith et al. 2005). Thus a 
model based only on the analysis of arc colors and their possible 
associations can lead to results that are not robust. A more de- 
tailed study of the color in the arcs, divided into components, is 
not the solution and may be wrong because the errors obtained 
from this kind of measurements beco me greater as we divide 
the arc into small fractions. For instance. lAlardl (12009) separated 
the individual elements in the arcs into two kinds: areas with 
bright sharp details and areas with smoother images. He found 
that bright regions have the same colors in the three arcs (as well 
as the smoother areas), and that the differences in colors between 
both areas (i.e. comparing the brightest ones with the smoother 
ones) were statistically significant. Thus, he argued that one sin- 
gle source produces the arc configuration in SL2S J02 140-0532. 
However, the magnitudes obtained by performing photometry in 
regions of 5 pixels in radius could have errors of ~0. 1 mag in 
bright details, and reach ~0.3 mag in the smoother areas where 
the signal is very low. Even with the ACS Snapshots, the errors 
are substantial, the images are not deep enough for this kind of 
analysis, and we need to be very cautious about his conclusions. 
As we discussed above, a model based only on colors would not 



be robust. This system clearly shows the need for a spectroscopic 
redshift to obtain an accurate lens model. 



2.3. Photometric redshifts 

Using the magnitudes reported in TableQ] we estimated the pho- 
tometric redshifts for the feat ures labeled A, B, and C with the 
HyperZ software (Bolzonella et al. 2000). We present the output 
probability distribution function (PDF) from HyperZ in Fig. [3] 
Arc C is constrained to be at z p h t = 1.0 ± 0.06, which is in very 
good agreement with the identification of the emission line as 
[OII]/B727 at z spec = 1.023 + 0.001. The multiply imaged sys- 
tem represented by arcs A and B has z p hot= 1.7 ± 0.1 and z p hot= 
1.7 ± 0.1, for the arcs respectively. Although the probability dis- 
tributions in the three cases are very broad (which is common in 
photometric redshift estimations using only five filters), there is 
no overlap in the solutions between system AB and arc C below 
a 3<x value. Along with the similarity between the distributions 
of arcs A and B, and the spectroscopic data, the possibility that 
A, B, and C belong to the same source is again ruled out. 

2.4. Pixelizing the arcs: increasing the number of constraints 

When a source is mapped onto a system that shows multiple sub- 
components (surface brightness peaks), similar to the one dis- 
played by SL2SJ02 140-0535 (system AB), the number of con- 
straints can be increased as well as the degrees of freedom (for a 
fixed number of free parameters), since the peaks can be conju- 
gated as different multiple image systems. To highlight the sub- 
structures appearing in system AB, we subtracted the six central 
galaxies of the group (Gl, G2, G3, and the three small ones near 
those, see Figure[T]top left). Following McLeod et al. ( 1998), we 
analyzed the F475W image and fitted a model convolved with a 
PSF We used six de Vaucouleurs profiles to fit the galaxies and 
a synthetic PSF because there was no suitable non-saturated star 
in the field of view. To account for possible differences, we also 
convolved the PSF with a Gaussian profile. In Fig. [4] we show 
the central region of the group after the galaxy subtraction. 

To avoid introducing any bias into the modelling, we conju- 
gated any bright peaks for which we were high confidence. This 
procedure transforms our simple AB system into four different 
systems, l.i, 2.i, 3./, and 4./, where i goes from 1 to 3. In Fig. [4] 
we show the image identification of the peaks, and in Table [2] 
we report the positions (we assume that the typical uncertain- 
ties in the positions of identified conjugated images are 1 pixel 
~ 0.05")- To keep the same notation and avoid confusion, we re- 
named arc C as 5. 1 . Therefore, our model have five different arc 
systems in the optimization procedure. 

In summary, we have shown that features presented in sys- 
tem AB originate from the same background source, which is 
multiply imaged. This background source at z p h t = 1.7 ± 0.1 
yields three images in the image plane: arc A composed of two 
merging images and arc B. For each image, we have been able 
to conjugate substructures (systems l.i, l.i, 3.;, and AS). In addi- 
tion, we have convincingly shown that arc C (now system 5. 1) is 
a different object that is not multiple imaged. 



3. Dynamical constraints 

We describe the profiles used in our strong lensing model and 
also present how the dynamical constraints enter into the model 
on both the galaxy and group scales. In short, we propose that the 
dark matter component of the group consists of a single large- 
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scale clump and we add smaller-scale clumps as galaxy-scale 
perturbations that are associated with the three individual galax- 
ies at the center of SL2S J02 140-0535. For these smaller clumps, 
we use the internal velocity dispersion of galaxy Gl as well as 
common scaling laws that links the parameters of the galaxies 
to their luminosity. In the large scale-clump, the dynamical in- 
formation obtained from spectroscopic data is used to set the 
possible range of values in the parameters that characterize the 
profile. 

3.1. Galaxy scale 

We assume that the density profile of the stellar mass distri- 
bution in the central galaxies of SL2S J02140-0535 follows a 
pseudo isothermal elliptical mass distribution (PIEMD). This 
profile was derived by Kassiola & Kovner ( 1993) and is param- 
eterized by a central density po, linked to the central velocity 
dispersion, <x , which in turn is related to the depth of the po- 
tential well (see lLimousin et al1|20 05; Eliasdottir et al. 2003, for 
a detailed discussion of the properties of this mass profile). It is 
described using two characteristic radii that relate to changes in 
the slope of the density profile. In the inner region, the profile be- 
haves like a non-singular isothermal profile with central density 
po and core radius r core . In the outer parts (r cl „ < r), the density 
progressively falls from pa r 2 topoc r~ 4 , introducing a natu- 
ral cut-off. A clump modeled using PIEMD can be characterized 
using seven parameters: the center position, (X, Y), the ellipticity 
e, the position angle 8, and the parameters of the density profile, 
Co* f core , and r cut . Below, we discuss the parameters associated 
with the dynamics of the central galaxies, namely, <xo and r cut 
and we return to those remaining in the next section. 

Our spectroscopic data include the two brightest central 
galaxies in the group, Gl and G2. However, given the poor 
quality of the G2 galaxy spectra (low S/N ratio), we only mea- 
sured the line-of-sight internal velocity dispersion of galaxy 
Gl. Using the penalized pi xel-fitting method developed by 
ICappellari & Emselleml (12004 . we found that cr* os =215 ±34 
km s" 1 for the G-band absorption line profile. Hereafter we de- 
note with an asterisk quantities associated with galaxy Gl. This 
velocity corresp onds to cr^ - 253 km s -1 , according to the rela- 
tion reported bv lEliasdottir et all (|2007) who studied the internal 
velocity dispersion of some galaxies in the Abell2218 cluster. 
Assuming that the galaxy profiles can be described by spheri- 
cally symmetric PIEMD profiles, Eliasdotti r et al.1 (120071) found 
that the ratio of the measured velocity dispersion (which is de- 
rived using spectroscopic data) to the fiducial velocity dispersion 
of the PIEMD profiles, (i.e., the ratio of the stellar to dynamical 
mass) has an almost constant value (<x* « erf /0. 85) across the 
region where the spectra were obtained. 

Since the total stell ar mass of the galaxy is related to <x* and 
r* ut by the expression (Eliasdot tir et al.ll2007l) 



M. 



1.57T 



' euro 



(1) 



the parameter r* ut can be constrained if we determine M, in an 
independent way. To measure the stellar mass of galaxy Gl, 
we followed two different methods. First, we converted the g- 
band magnitude into an absolute rest frame V-band luminosity 
and found that Ly - 3.0xl0 10 L© v . We then used the mass-to- 
light ratio reported in SLACS lenses dGavazzi et alJ 12007 ) to 
estimate the stellar mass and found that M* ~ l.OxlO 11 M©. 
For the second metho d, we employed the HyperZ software 
(Bolzonel la et alJ l2000) with the five magnitudes of galaxy Gl 
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Fig. 4. F475W image ( 22"x22") after subtraction of the six cen- 
tral galaxies of the group. The green circles show the positions 
of the conjugated images that constitutes the AB arc system. The 
brown circle, mark the location of arc C. 



Table 2. Multiple image identification. 



Systems 


(!' 


6 




(J2000) 


(J2000) 


1.1 


02:14:07.684 


-05:35:27.37 


1.2 


02:14:08.184 


-05:35:27.54 


1.3 


02:14:08.477 


-05:35:31.68 


2.1 


02:14:07.708 


-05:35:27.19 


2.2 


02:14:08.165 


-05:35:27.32 


2.3 


02:14:08.493 


-05:35:31.69 


3.1 


02:14:07.739 


-05:35:27.26 


3.2 


02:14:08.110 


-05:35:26.95 


3.3 


02:14:08.517 


-05:35:32.00 


4.1 


02:14:07.770 


-05:35:27.35 


4.2 


02:14:08.048 


-05:35:26.93 


4.3 


02:14:08.538 


-05:35:32.48 


5.1 


02:14:08.015 


-05:35:38.54 



reported in Table Q] We found a slightly greater value of M* w 
4.0x10" M©. Once we apply Eq.Q] these range of masses pro- 
duce a range of possible values for r* cut of 1 - 6 kpc. 

The parameters of the two remaining galaxies are obtained 
using the scaling relations 

r , = r* (-Ml/2 

'cut — • cut\L"' ' 

cro = <r*(i0 1/4 , 

where cr^, r* ut and L* are set by those values of galaxy Gl dis- 
cusse d above. For a discussi on of these scaling relations, we re- 
fer to Limousin et all d2007l) . 



(2) 
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3.2. Group scale 

3.2.1. Velocity dispersion and virial mass 

In Fig. [5] we show the redshift distribution of the galaxies in the 
field of SL2S J02140-0535; as we can note, the spectroscopy re- 
veals a well defined group without any evidence of bimodality. 
The dynamical study of this group will be presented in a forth- 
coming publication (Munoz et al. 2010 in preparation) as part 
of a detailed study of seven strong-lensing galaxy groups in the 
SL2S. Here we briefly describe the method used. We adopted 
the formalism by Wilman et al. (2005) in order to determine the 
group membership. We identified the group members as follows: 
the group was initially assumed to be located at the redshift of 
the main bright lens galaxy, zi ens , with an initial observed-frame 
velocity dispersion of cr ;, s = 500(l+z/ e „. s )kms _1 . After comput- 
ing the required redshift r ange for group membership and apply- 
ing a biweight estimator (Be ers et aLl ll990) the iterative process 
was found to have reached a stable membership solution with 16 
secure members and cr = 630 ± 107 kms -1 . These galaxies are 
shown with red rectangles at the bottom of Fig. Q] 

Assuming that the group has an isotropic velocity dispersion 
and is no longer undergoing net expansion or contraction, we 
apply the virial theorem to estimate the mass of SL2S J02140- 
0535. Since in our lensing models we probe the two-dimensional 
(2D) mass of the groups, we define M v - 3nR v cr 2 o J2G as the 
virial mass inside the projected radius R v . Using the measured 
crios - cr/ V3 = 364 kms _I and a virial radius R v - Rh — 0.8 ± 
0.3 Mpc, where R h is the projected harmonic mean radius (e.g. 
Ilrgens et alj|2002l) . we obtain a mass of (1.1 ± 0.8) x 10 14 M©. 
We stress that this is not the true 2D mass at R v , but only an 
estimate of the mas s inside R v , and th at this mass is probably 
underestimated (e.g. lBiviano et alj|2006l) . 

3.2.2. Dynamical constraints in the NFW profile 

We describe our gal axy-group density profile as an NFW profile 
dNavarro et alJll997h 



P(r) 



Ps 



(r/r s )(l+r/r s ) 2 



(3) 



where p s is a characteristic density, and r s is a scale radius that 
corresponds to the region where the logarithmic slope of the den- 
sity equals the isothermal value. To link the velocity dispersion 
calculated in the last subsection with the velocity associated with 
the potential of the NFW halo, we follow two different and in- 
dependent methods. First we define the velocity in terms of the 
scale radius and the characteristic density to be 



cr^^l+luiWp, 



(4) 



It is defined in such a way that it represents a realistic veloc- 
ity dispersion and not only a scaling parameter, which is useful 
when comparing the velocity predicted by strong lensing models 
with the velocity dispersion measured in dynamical studies (see 
the discussion in Verdugo et al. 2007). In a second way, follow- 
ing lEliasdottir et alj (J2007I) we calculate the projected average 
line-of-sight velocity dispersion inside the radius R v as 



, 2n $ * oi'K)Y,'R)RdR 



«.X*v) 



M 2D (R,) 



y 












1 






■ 






SLi!SJUai4UB-UMt>32 






; 









z = 0.445 




^ 4 




A 




- 


ti 


0.43 


0.44 0.45 


0.46 


6 










- 


4 










- 


2 


- 


1 


| 










- 


n 




. 


1 


.1 11 


w 


, 1 


, 1, 1 





0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



Fig. 5. The histogram (bin size A = 0.004) shows the redshift 
distribution along the direction of SL2S J02140-0535. The spec- 
troscopic data reveals a well-defined group at z — 0.444. The 
inset highlights the central part rebinned to A = 0.0025. There 
are 17 galaxies in the histogram, but after the iterative algorithm 
(see text) one outlier was removed. 



where R is a 2D radius and cr 2 (R) is the projected line-of-sight 
velocity dispersion given by 



(6) 



2 2G r^uifW)^—^, 

and X(/?) the surface mass density. We note that Eq. [5] depends 
not only on the virial radius, R v , but also in the parameters that 
characterize the profile. Integrating Eq. [3] we obtain 



M m (r) = 4«7^p a 



ln(l+r/r,)- 



r/r s 



1 + r/r s 



(7) 



Integrating the surface mass density (see Appendix), we then 
obtain the expression for M2d(R\), and substituting Eqs.|7]and|3] 
into Eq.|6] it is possible rewrite Eq.[5]as 



G(R v ,r s ) 



(crl s ){R v ) = 2nGp s ^ 



F(R v ,r s ) 



(8) 



where 



G(R 



,r,) = 

JO JR 



ln(l +r/r,)- 



r t + r 



1 



Vr 2 - R 2 



(r/r.Xl+r/r,) 2 



-RdrdR 



(9) 



and F(R V , r s ) is given by Eq. IA.4l We emphasize that the depen- 
dence of the projected average line-of-sight velocity dispersion 
(Eq. [8} on the virial radius, the scale radius, and the characteris- 
tic density. If we assume that Eqs.|4]and|8]are related by the ex- 
pression (erf )(R V ) = cr 2 /3 (which is qualitatively correct in the 

sense that cr represents a 3D velocity dispersion and (erf )(R V ) 
represents the projected average line-of-sight velocity dispersion 
assuming no anisotropy and a spherically symmetric NFW pro- 
file), then 



(5) 1 



3n 



G(R v ,r s ) 



2(l+]nl/2y sF(Rv ,r s ) 



(10) 
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Thus, given R v , we can use the above equation to obtain the pos- 
sible range of values for the scale radius. It is noteworthy that 
although the velocity dispersion does not appear explicitly in 
the equation, the equation is not independent of this value since 
our virial radius is computed from the relative positions of the 
confirmed members, which in turns depends on the velocity dis- 
persion of the group. In other words, the velocity dispersion is 
implicit in the equation, and therefore, in our calculations of the 
scale radius. For R v - 0.8 ± 0.3 Mpc we obtain r s - 150 ± 50 
kpc. It then follows that c^fw - Rvfrs =5 + 3. After the scale 
radius is estimated, the 2D mass of the halo expressed in terms 
of the velocity dispersion of the group and the virial radius can 
be determined in a simple fashion using Eqs.l4land lA.3l and ex- 
pressed as 



M 2D (R V , r s , o-) 



2na r s 



(1 +lnl/2)G 



F(R v ,r s ). 



(ID 



Using our previously computed values, this yields M 2 d = (4 ± 
2) x 10 14 M0. This estimate is, within the errors, slightly greater 
than the value M v computed from dynamics, although we show 
in Section [5] that it is in good agreement with the strong and 
weak lensing mass estimates. For the strong lensing analysis (see 
Section |4), we use these dynamical results to put some priors in 
the range in which the parameters are allowed to vary: we set the 
velocity dispersion range to be 523 km s < <x< 737 km s -1 , and 
the scale radius to be 100 kpc < r s < 200 kpc, i.e. to vary within 
1-cr of their values. 



4. Gravitational lensing analysis 

4.1. Strong lensing 

To study the lensing group SL2S J02 140-0535, we used the para- 
metric rnethod_jmplemented in the LENSTOOLQ ray-tracing 
code (Kn eiblll993l) . This software use a Bayesian Monte Carlo 
Markov chain (MCMC) method t o search for th e most likely pa- 
rameters in the lens modeling (Juli o et al.ll2007l) . 

The identification of the multiply imaged system presented 
in Section|2] leads to 16 observational constraints. With these 
constraints, we computed a model (optimized in the image 
plane) with the following set of free parameters {X, Y, e, 6, r s , 
<r, zi, zi, Z3, Z4, r* ut }, where the first six parameters characterize 
the NFW profile and Zi represents the redshifts for the systems 
\.i, 2d, 3,i, and A.i, respectively. All the parameters are allowed 
to vary with uniform priors except r s and <x for which we used 
Gaussian priors with the 1-cr errors computed in the last section. 
System 5.1, at z S pec - 1-023 ± 0.001 is used in the analysis as 
a singly imaged object. Multiple image systems were allowed 
to vary between 1.5 and 2.0 (approximately 2-<x from the value 
inferred from the photometric redshift analysis) to take into ac- 
count the broad PDF common to these measurements; and for 
the same reason we do not choose a Gaussian prior in this case. 
The remaining parameters, those that describe the central galax- 
ies, were set as follows: the center of the profiles and the ellip- 
ticity and position angle are assumed to be the same as for the 
luminous components. We use the values reported in Table Q] 
The velocities dispersions are given by Eq. [2] for a <x* = 253 

kms -1 . 

Figure [6] shows the PDF for the arc systems, and the best- fit 
model solutions for z, obtained after the optimization. The mean 



This 



software 



http://www.oamp.fr/cosmology/ lenstool/ 
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Fig. 6. Normalized PDFs. The red, blue, green and magenta his- 
tograms shows the distributions obtained with our model, for z\, 
Zz, Z3 and z\ respectively. The vertical lines corresponds to the 
best-fit model solutions given by the optimization procedure. 
Black continuous and black dashed lines are the distributions 
shown in Figure [3] for arcs A and B, respectively. 
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Fig. 7. Tangential critical lines (red) and caustic lines (blue) for 
two different sources located at z — 1.0 and z - 1.7. As in Fig. [4] 
the circles show the measured positions of the images (input data 
for the model). The squares (orange and cyan) are the model- 
predicted image positions after optimization, and the asterisks 
(orange and cyan) the positions in the source plane. Although the 
cyan asterisk appears to over lay the caustic line, this is only an 
artifact, a product of the thick lines, and the size of the asterisks 
that we used to highlight the figure. The source is not over the 
caustic line. 



optimal value is z - 1.7 ± 0.2. It is clear from the figure that the 
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Table 3. Best-fit model parameters. 



Comp. 



X 



Y 
["J 



[°] 



[kpc] [kpc] 



[kpc] 



c 200 CO °~ 

[kms -1 ] [kms 4 ] 



Xdof 



1.2 



Group 1.3 ±0.1 0.5±0.3 0.28 ± 0.02 111.6 + 0.2 — 170 ± 18 — 6.0±0.6 — 664±18 

L' — [0.15] — 5.6 ±1.3 [253] — 



Notes. The first column identifies the different scale components. The L* denotes the galaxy-scale mass component. Columns 2 and 3 list the 
position in arcseconds relative to the BGG. In columns 4 and 5 we provide the geometric parameters. From columns 6 to 11, we present the 
different profiles parameters, and in the last column Xhe-Xoop- Values in square brackets are not optimized. 



model PDFs are slightly far to be Gaussians and z values greater 
than two are still favored. This reflects the lack of arc systems 
(only one system with a fixed redshift) available to constrain the 
model and shows the importance of using the photometric red- 
shift as a prior. The results of our best fits are summarized in 
Table[3] In Figure|7] we show the predicted positions of our best 
model (orange and cyan squares), as well as the observed posi- 
tions (green and brown circles). The remarkable agreement be- 
tween both positions is quantified by a mean scatter in the image 
plane smaller than 0.04". We obtained Xdof ~ ^' snowm g that 
the extent of the match between the observed positions and those 
predicted by the model is in accord wit h the error varian ce. This 
X^of * s defined in the image plane (see lJullo et al.ll2007l) assum- 
ing the same uncertainties in the positions for all the images (see 
Sect. l2.4l >. We note that our model do not predict extra images, 
which supports the reliability of our strong lensing model. 



4.2. Weak lensing 



T he weak lensing analys is of this galaxy group was performed 
in iLimousin et al.l ([2009). We refer the interested reader to that 
paper for a detailed description of the methodology, and present 
below a brief summary. They selected as background sources 
the galaxies whose i band magnitude falls between 21.5 and 
24. Their density equals to 13 arcmin~ 2 . The completeness 
magnitude in this band is 23.96 and the seeing equals 0.61". 
A Bayesian meth od, implemented in the Im2shape software 
(Brid le et al.l2002l) . is used to fit the shape parameters of the faint 
background galaxies and correct for PSF smearing. From the cat- 
alogue of background galaxies, Limo usin et al.1 (2009) pursued 
a one-dimensional weak lensing analysis. They fitted a singular 
isothermal sphere model (SIS) to the reduced shear signal be- 
tween 150 kpc and 1.2 Mpc from the group center, finding an 
Einstein radius equals to 3.6 + 2.4". To relate the strength of 
the weak lensing signal to a phy sical veloc ity dispersion char- 
acterising the group potential, ILimousin et al.l d2009) estimated 
the mean geometrical factor using the photometric redshift cat- 
alogue from the T00 04 release of the CFHTLS-Deep survejy 
dlenna & Pell6ll2006l) . They found that cr SIS = 612+^° kms 1 . 
This translates into a projected mass computed within a circular 
aperture of radius equals to 1 Mpc of (2.8+1.5) x 10 14 Mq. 



http://www.ast.obs-mip.fr/users/roser/CFHTLS_T0004/ 



5. Discussion 

Our best-fit model demonstrates that it is possible to re- 
cover the observational features presented in the galaxy group 
SL2SJ02 140-0535 using a single halo for the group compo- 
nent plus three galaxy scale components associated with the cen- 
tral galaxies. We do not need extra assumptions about the mass 
distribution in the group. In addition, the position angle of the 
halo is oriented with a well-constrained position angle at 9 - 
1 1 1.6+0.2, which is the same direction as the one defined by 
the luminosity contours (Limousin et al. 2009) and is consistent 
with the spatial distribution of the confirmed members of the 
group (see Fig. Q] bottom). This supports a scenario where the 
mass is traced by light. 

The degeneracies in the <T-r s space are very common in 
strong lensing modeling (e.g. lJullo et al.ll2007T) since the region 
where the arcs appear is limited to the central region of the clus- 
ters and the scale radius is generally three or four times the dis- 
tance to that region. In galaxy groups, this region is indeed, even 
smaller than in clusters. Thus, it is impossible to constrain the 
scale radius of a group using only its observed strong lensing 
features. To highlight this behavior, we depicted in the top panel 
of Figure [8] the PDFs for the cr and r s parameters. We note how 
models with a scale radius greater than 200 kpc are still possi- 
ble (in addition to those smaller than 100 kpc), but we placed a 
limit on this value using information from the group dynamics. 
Thus, this give us a method that would be useful in breaking the 
degeneracies in the cr-r s space using the dynamics of the group 
obtained from spectroscopic data. This is important because it 
provides a technique to fill the gap between strong lensing and 
weak lensing constraints in galaxy groups. 

To analyze the dark matter profile of SL2SJ02 140-0535 in 
greater detail, we calculated the concentration C200 from the val- 
ues of r s and cr (see the discussion in Verdugo et al. 2007). Here 
we assumed that C200 = r 200 lr s , where we defined r2oo as the ra- 
dius of a spherical volume within which the mean density is 200 
times the critical density of the universe at the given redshift z 
of the group. The bottom panel of Figure [8] shows the 2D his- 
togram for the parameters C200 and r s , where we can note the 
same tendency with r s as in the top panel. Our best-fit model 
infers that C200 = 6.0 ± 0.6, which is in good agreement with 
the value Cnfw presented in Section [4] On the other hand, if we 
compare it with the values expected from ACDM numerical sim- 
ulations, w e found that our con centration is slightly higher. From 
the work of iDuffv et all (120081) . a DM halo at z = 0.44 with M 200 
* 1 x 10 14 M has C200 ~ 4.0. 
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Fig. 9. The 2D projected mass as a function of the aperture ra- 
dius measured from the BGG. The green and blue shaded areas 
corresponds to the mass profile within \cr errors for the strong 
lensing and the weak lensing models, respectively. The red re- 
gion illustrates the 2D mass derived from Eq. QT| as a function 
of R v , The point with error bars is the mass estimate from dy- 
namics. The orange-shaded region shows the area where the arc 
systems lie and the gray region shows the model-predicted scale 
radius. We also depict with an arrow the zone in which the weak 
lensing signal was measured. 



nal strong enough to allow us to fully p robe an NFW profile. 
Until the work of Limousin et al. (2009), weak lensing signals 
on these mass scales had indeed been reco vered only by stack- 
ing group images (e.g. Mandelbaum et al. 2006). However, the 
strength of the signal in SL2S J02 140-0535 is sufficient to fit a 
SIS and obtain the projected mass as a function of the radius. 
In Figure [9] we compare the masses obtained with the differ- 
ent methods adopted in the present work: strong lensing (green), 
weak lensing (blue), M2D derived from Eq. QT| (red), and the 
virial mass (point with error bars). As we place limits on the 
possible values of the scale radius, the strong lensing mass is still 
reliable up to the scale radius (r s - 170 ± 18 kpc), but beyond 
this value the extrapolation of the mass becomes less secure. 
However, from Figure [9] we can note that the first three masses 
are in agreement within the errors in the range 0.5 Mpc < r < 1.1 
Mpc, i.e. the interval defined by R v . The virial mass, despite be- 
ing lower than all masses, differs more from Mid- However, we 
stress again that it is not the true 2D mass at /?,, but only an esti- 
mate of the mass inside R v . This mismatch between virial mass 
and lensing mass is also found in other groups, but the analysis 
of this subject is beyond the scope of the present work and will 
be discussed in detail in a forthcoming paper (Mufioz et al. 2010 

in preparation). 

Using a perturbative method, Alarg (120091) reconstructed the 
structure and potential of the gravitational lens SL2S J02140- 
0535. In his work, he assumed that the A, B, and C arcs belong 
to the same source. He argued that the differences in colors ob- 
served in the HST-ACS images were related to the crossing of 
the caustic line by part of the source. His model predicts that 
one half of image C comes from the same part of the source as 
images A and B, but the other half corresponds to an area of the 
source that has no counterpart in the large arcs (A and B). The 
perturbative method is capable of reproducing this arc configu- 
ration, but the model predicts the existence of an independent 
dark component that does not follow light. In the present work, 
we have demonstrated that AB system and C system belong to 
different sources, one at z — 1.7 ± 0.1 and the other one at z S pec 
= 1.023 ± 0.001. Our model, with a single NFW profile (which 
is consistent with the analysis of the luminosity contours, Fig.Q] 
bottom, and the redshift distribution of the galaxies in the field 
of SL2S J02 140-0535, see Fig. EJ, provide a good fit of the arc 
features observed in the group and our results point out that the 
projected light of the group traces the dark matter distribution. 



6. Conclusions 

We have probed the gravitational potential of SL2SJ02140- 
0535 on a large radial range using complementary techniques: 
strong lensing, weak lensing, and dynamics. In the strong lens- 
ing regime, we performed a fit using an NFW profile and three 
galaxy-scale mass components as perturbations to the group po- 
tential. We used PIEMDs for the individual galaxies and set 
constraints on their geometrical and dynamical parameters us- 
ing observational data (image and spectroscopy). We measured 
the velocity dispersion of the group using Keck (LRIS) and VLT 
(FORS2) spectra, and we applied this information to constrain 
and build a reliable strong lensing model based on one multiply 
imaged system and one single image syst em. In the weak lensin g 
regime, we used the analysis reported in lLimousin et al.l (|2009). 
Our main results can be summarized as follows: 



Given the relatively low mass of galaxy groups (compared 
with clusters), we do not expect to derive a weak lensing sig- 



1 . We have shown that the AB and C systems belong to differ- 
ent sources, one at z — 1 .7 + 0. 1 (photometric redshift) and 
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the other one at z — 1.023 ± 0.001 (spectroscopic), respec- 
tively. 

2. Our best-fit lens model reproduce quite well the image sys- 
tems observed in the field of SL2S J02 140-0535, with no ex- 
tra assumptions about the underlying mass distribution (e.g. 
bimodal distribution). This agrees with the spectroscopic 
analysis that shows a well-defined group. Furthermore, we 
obtained a position angle of the halo of 6 — 111.6+ 0.2, that 
is consistent with the orientation of the luminosity contours 
and the spatial distribution of confirmed members. 

3. The spectroscopic information provide constraints on the 
scale radius of the NFW profile. In the present work, we 
showed that is possible to use the dynamics at group scale 
to discriminate between a great range of possible r s values 
and cr (in particular 100 kpc < r s < 200 kpc and 523 km s _1 
< cr < 737 kms _I , respectively). 

4. The scale radius and the velocity dispersion found in our 
best-fit models, r s - 170 ± 18 kpc and cr - 664 ± 18 km s -1 , 
respectively, yields a concentration value C200 = 6.0 + 0.6, 
which is consistent with c^fw — 5 ± 3, but slightly greater 
that the one predicted by ACDM simulations. 

5. The masses obtained with the strong lensing, weak lensing, 
and Mid methods agree within the errors, but the virial mass 
is slightly below these values. Thus, we demonstrated that it 
is possible to combine these different methodologies to get a 
complete insight into the mass of a galaxy group. 



where 



F{® = 



arccoshj I if £ < 1 




(A.2) 



and the dimensionless coordinate £, is the radius in the XY plane 
in units of the scale radius, £, - (x/r s , y/r s ). It then follows that 
the 2D mass inside the virial radius can be expressed as: 



M 2D (p s ,R v ,r s ) = inp s ijF(R v ,r s ), 
where 



(A.3) 



F(R v ,r 



Wo JF^i 



-^arccoshl 

Vw 1 * 



r* 1 1 1 

I — ^=^ 1 ^=^= arccos - 



d%. (A.4) 



and £ = R v /r s . We note that the mass depends on both the char- 
acteristic density and the scale radius of the NFW profile. 



The lensing and dynamical analysis of SL2SJ02 140-0535 
presented in this paper, has had two principal points. First, it has 
investigated the impact of including spectroscopic information 
of arcs (i.e. secure arc identifications) in strong lensing models 
to avoid complex mass distribution models that do not follow the 
light. Secondly, we studied the possibility of constraining large 
scale properties that could not be measured by strong lensing, 
such as the scale radius in the NFW profile. Future spectroscopic 
follow-up of the arcs A and B will provide confirmation of our 
model. In addition, the deepest HST images of the group will 
allow us to construct more accurate weak lensing models capable 
of probing an NFW profile on large scales, as well as improving 
our strong lensing model. 
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Appendix A: M 2 /> inside the virial radius R v in the 
NFW profile 

The surface mass de nsity in the NFW profile is given by 
(iGolse & Kneibll2002l) 



liO = 2p s r s F(&, 



(A.l) 
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